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Effects of orbital degeneracy on magnetic susceptibility in paramagnetic phases are inves- 
tigated within a mean-field theory. Under certain crystalline electric fields, the magnetic mo- 
ment consists of two independent moments, e.g., spin and orbital moments. In such a case, 
the magnetic susceptibility is given by the sum of two different Curie- Weiss relations, leading 
to deviation from the Curie- Weiss law. Such behavior may be observed in d- and /-electron 
systems with t2 g and F$ ground states, respectively. As a potential application of our theory, 
we attempt to explain the difference in the temperature dependence of magnetic susceptibilities 
of UO2 and Np0 2 . 
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The Curie- Weiss law is a well-known fundamental re- 
lation that characterizes the magnetic response of con- 
densed matter. 1 Indeed, in various kinds of materials, 
magnetic susceptibility \ is known to follow the Curie- 
Weiss law given by 

X = C/(T-9), (1) 

where T is the temperature, C is the Curie constant, and 
6 is the Weiss temperature. The Curie- Weiss law is ob- 
tained by applying mean-field theory to a model in which 
magnetic moments interact with one another under an 
applied magnetic field. In the mean-field approximation, 
C oc /Lt ef j and 9 — CA, where /z e ff is the effective magnetic 
moment and A is the coupling constant of the interaction 
between magnetic moments. Thus, /x e fi and A can be es- 
timated by fitting the experimental results for x to the 
Curie- Weiss form. 

On the other hand, there exist many materials in which 
the magnetic susceptibility simply does not follow the 
Curie- Weiss law. One reason for this, of course, is the ef- 
fects of fluctuations, which are not included in the mean- 
field approximation. In order to improve the theory, one 
must take into account the effects of fluctuations beyond 
the mean-field approximation in the temperature region 
near a phase transition point. 

However, even if the temperature is so high that fluc- 
tuation effects can be ignored, the deviation from the 
Curie- Weiss law is still significant in some materials with 
active orbital degrees of freedom. In such a case, the de- 
generacy of the ground-state multiplct is generally lifted 
by the effects of a crystalline electric field (CEF) . When 
the interactions between magnetic moments are consid- 
ered within a mean-field theory, we obtain the magnetic 
susceptibility as 

X' 1 = Xcef - ^ ( 2 ) 
where xcef indicates the magnetic susceptibility of a sin- 
gle ion in a CEF. When the CEF level splitting between 
the ground and excited states is much larger than ther- 
mal energy, it is enough to consider only the CEF ground 
state to evaluate xcef- In this case, since xcef follows 
the Curie law, that is, xcef = C/T, \ again follows the 
Curie- Weiss law. 



From the above discussion, when the temperature is 
higher than the phase transition point and the excited 
state is well separated from the ground state, a Curie- 
Weiss law seems to hold, but it is premature to conclude 
it. An important point has been missed. Thus far, we 
have implicitly assumed a one-component magnetic mo- 
ment. In some materials with an active orbital degree 
of freedom, however, the magnetic moment consists of 
two components. For instance, readers can envisage a d- 
electron system in which orbital angular momentum is 
not quenched. In such a system, the magnetic moment 
is given by M = —/.ib(2S + L), where /j,b is the Bohr 
magneton, S is the spin angular momentum, and L is 
the orbital angular momentum. 

In this Letter, we study magnetic susceptibility in the 
paramagnetic phase of a cubic system. Note that in such 
a system with high symmetry, the orbital degree of free- 
dom becomes active in general, and the magnetic mo- 
ment consists of two independent moments in some cases. 
It is intuitively understood that in such a situation, the 
magnetic susceptibility is given by the sum of two dif- 
ferent Curie- Weiss relations, leading to non-Curie- Weiss 
behavior. This is our main message, but an explicit form 
of the magnetic susceptibility of the system with two 
independent moments is needed for the analysis of ac- 
tual materials. Thus, first we derive magnetic suscepti- 
bility from a phenomenological model in a level of stan- 
dard mean-field approximation by using simple algebra. 
As examples, we consider a d-electron system with a t2 g 
ground state and an /-electron system with a L 8 CEF 
ground state. In particular, we compare our magnetic 
susceptibility with the experimental results of Np02 2 ~ 4 
with a T 8 CEF ground state 5 ' 6 and U0 2 7 ' 8 with a T 5 
CEF ground state. 9, 10 The difference in the temperature 
dependence of magnetic susceptibilities in these materi- 
als is naturally explained by considering the fact that 
dipole and octupole moments coexist in Np02, whereas 
only the dipole moment exists in UO2. 

To consider both situations with and without spin- 
orbit interaction, we generally express the magnetic mo- 
ment M as M = M (1) + M (2) in the CEF ground- 
state multiplet, where the operators and 
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are assumed to be orthogonal to each other, namely 
tr (Ma ' 1 Afi 2 ') = 0. Here, a indicates a cartesian compo- 
nent and tr(- • • ) denotes the sum of expectation values 
over the CEF ground states. For simplicity, we further 
assume that a magnetic field H is applied along the z- 
direction. The effective Hamiltonian, then, is given by 

H=-HY J M iz - £ lt%M^M$, (3) 

where i and j denote lattice sites, and I^jp is the inter- 
action between M-^ and M^J . Note that the Hamil- 
tonian should also include more complicated interac- 
tion terms among the multipolc moments, but in the 
mean-field approximation for a paramagnetic phase up 
to 0(H 1 ) 1 such terms vanish. Thus, they are simply 
dropped at this stage. 

Let us now apply the mean- field approximation to TL. 
Since we are considering the magnetic susceptibility in 
the paramagnetic phase, the site dependence can be sup- 
pressed, that is, = M„ M ' . In addition, for a / ft 
we can show J2j lj£ "jp = by using a symmetry argu- 
ment about a cubic lattice. For example, when the vector 
connecting sites i and j± is given by r± — (±x,y,z), 
the relation = holds in a cubic lat- 
tice, indicating that J2j ^ix 'jy = 0- Namely, the mean- 
field from to should be zero. As a result, 

(M^) = (M^ ] ) = under a magnetic field H along the 
z-direction, where (• • • ) denotes the expectation value. 
Thus, the mean-field Hamiltonian for the paramagnetic 
state up to 0(H 1 ) is expressed by 

Hmf = -HMt-^MPiMP) + X 2 MP(M^)}, (4) 

where A„ = 2 J2i j Ii z $ z /N an d N is the number of sites. 
Note that we have chosen the decomposition of the mag- 
netic moment so as to satisfy A12 = 2 J2i j H z j z /N = ®- 
Note also that in general, there are other moments which 
couple with magnetic moments, 11 but symmetry argu- 
ments tell us that the mean-field from such moments 
cancels out in the paramagnetic state, as the cancella- 
tion of the mean-field from M^' to for a /ft On 

the other hand, M z ^ and M z 2 ^ have the same symmetry. 
Thus, eq. (4) does not change even if we include other 
interactions, unless we consider spontaneous symmetry 
breaking. 

The magnetic susceptibility for each moment, defined 
as Xv = d(M| !/ ^)/d_ff (y = 1 or 2), follows a Curie- Weiss 
law independently as Xi/ = C V /(T — 6 V ), where 9 V — 
\ V C V . Note that the Curie constant C v is not determined 
here, since we do not set the explicit form for M^K 
Then, the magnetic susceptibility is given by \ — X1+X2, 
but it is useful to express it as 

= T + 6-9,-62 

X (T-0 1 )(T-0 2 y 

where C = C\ + C 2 and 6 = (d6»i + C 2 2 )/C. At 
high temperatures, the magnetic susceptibility follows 
the Curie- Weiss law \ = C/(T — 9) asymptotically. We 
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Fig. 1. Temperature dependence of the inverse of the magnetic 
susceptibility \ (solid curve) obtained from eq. (5). The dashed 
line represents the Curie- Weiss behavior at high temperatures. 
The parameters are set as C\ = C 2 = C/2, 6»i = |0|, 8 2 = -3|0|, 
and < 0. 

note that 9 plays the role of the Weiss temperature for 
the high-temperature behavior. 

Figure 1 shows a typical example of the temperature 
dependence of x- First, we remark that x diverges at a 
temperature 9* , where 9* is defined as the larger of 9\ and 
9 2 , and it is always larger than or equal to 9. We also note 
that 9* can become positive, even if 9 is negative. This 
fact indicates that if we naively assume a Curie- Weiss 
law, we may obtain different Weiss temperatures. In the 
worst case, one may mistake even the sign of the Weiss 
temperature, which determines whether the interaction is 
ferromagnetic or antifcrromagnctic. Thus, we should pay 
due attention to this effect when analyzing the magnetic 
susceptibility of such systems. 

Let us now discuss how our theory works in actual 
systems. For this purpose, we choose two cases. One is a 
(i-electron system including only t 2g ground states, and 
the other is an /-electron system having a r 8 quartet 
ground state. 

First, we consider a t 2g system without the spin- 
orbit interaction. There is no linear coupling between 
the spin and orbital angular momenta due to the ro- 
tational symmetry in spin space. Thus, we can simply 
decompose M z = mP + ikfj 2) as M z 1] = -2fj, B S z and 
M z = —[j,bL z . In the t 2g subspace, the eigenvalues of 
the orbital angular momentum L z and the spin angular 
momentum S z are given by L z = —I, 0, 1 and S z = ±1/2, 
respectively, where A denotes the eigenvalue of opera- 
tor A. Then, we obtain eq. (5) with C\ = /u|/A;b and 
C 2 = 2 / u|/(3/cb), where fee is the Boltzmann constant. 

Note that we have considered an ideal model in order 
to see the role of the orbital angular momentum. To in- 
clude the effect of the spin-orbit interaction, we should 
change the decomposition of M z . Moreover, in actual 
materials, the t 2g level is easily split into sublevels by 
Jahn- Teller distortions, and the orbital angular momen- 
tum tends to be quenched. Thus, the observation of the 
effects of the orbital angular momentum on the magnetic 
susceptibility in t 2g systems may be difficult. Neverthe- 
less, efforts to realize such a possibility are believed to 
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Fig. 2. Schematic views for orthogonality between the dipole mo- 
ment J and Ta u octupole moment T. (a) In the Tg subspace, 
J and T are linearly independent, but not orthogonal to each 
other. Then, we define T' by subtracting the component parallel 
to J from T. (b) In the Ts subspace, T is parallel to J. 



be challenging and interesting. 

Next, we turn our attention to an /-electron system. 
As is well known, due to the strong spin-orbit interaction, 
the total angular momentum (dipole moment) J , rather 
than S or L, is a good quantum number in /-electron sys- 
tems. Thus, one is not allowed to decompose M into spin 
and orbital angular moments as in the i2g electron sys- 
tem. However, in /-electron systems, multipole moments 
higher than dipole can influence the magnetic suscepti- 
bility under certain CEF potentials. In fact, in a system 
on a cubic lattice with a Ts CEF ground state, the oc- 
tupole moment T with the same symmetry as the 
dipole moment J becomes a degree of freedom which is 
independent of the dipole moment. 12 Then, we can ap- 
ply our theory to such a dipole-octupole system. It is our 
message that the coexistence of dipole and octupole mo- 
ments can naturally lead to non-Curie- Weiss behavior in 
some /-electron systems. 

We describe an effective model for the dipole-octupole 
system. First, in the case of a strong spin-orbit inter- 
action, M is proportional to J through the well-known 
relation M — ~gj\x^J , where g,j is Lande's g-factor. 
Second, in the T$ subspace, we define an octupole mo- 
ment T' as a linear combination of the original r<i u oc- 
tupole moment T and J, in order that T' is orthogonal 
to J [see Fig. 2(a)], namely, tr(T^J a ) = 0. We normalize 
T' so as to satisfy tr(J Q J a ) = tr{T' a T' a ). In terms of J 
and T', the mean- field Hamiltonian is given by 

Hmf = - HM Z - (gj^ B ) 2 [X d J z {J z ) + KT' Z {T Z ) 

+ X do (J z (T' z )+T z (J z ))}, 



(6) 



where Ad, A , and Ado are coupling constants of dipole- 
dipole, octupole-octupole, and dipole-octupole interac- 
tions, respectively. 

Now, we define M z (1) and so as to diagonalizc 

eq. (6) into the form of eq. (4). Then, we can exploit 
the present result for magnetic susceptibility, indicating 
that Mz and follow the Curie- Weiss law indepen- 
dently. We obtain eq. (5) also for Tg quartet systems, 
with 0, 9i, and 02 satisfying the relations — CAd, 
0i + 02 = C(X d + A ), and X 2 = C 2 (X d X Q - X 2 do ). It 
has already been pointed out that the two moments 
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Fig. 3. Temperature dependence of the inverse of the magnetic 
susceptibilities \ of UO2 and of Np02- Solid curves indicate 
the experimental results for H = 10 kOe and H \\ [100] (a) 
U02- 8 The dotted curve represents the Curie-Weiss law with 
,u efI = 3.23mb [C = Mcff/( 3fc B)] and 9 = -223 K. (b) Np0 2 . 4 The 
dashed curve represents the fit obtained with eq. (5) with /i e ff = 
3.20mb, = -146.1 K, 6»i = -161.9 K, and 62 = -22.5 K. The 
dotted curve represents the Curie- Weiss fit with fi e ff = 3.20mb 
and 9 = -139 K. 



follow the Curie- Weiss law with different Curie and Weiss 
constants at least for certain interaction forms, 12 but 
here we emphasize that this property is retained in gen- 
eral when we choose the proper decomposition of M . 

As a typical application of eq. (5) to actual materials, 
let us discuss the magnetic susceptibility of NpC>2, 2 ~ 4 
in which the CEF ground state has been found to be 
a quartet. 5,6 For comparison, we also consider the 
magnetic susceptibility of UO2 7 ' 8 with a T 5 triplet CEF 
ground state. 9,10 Note that we can apply the present 
theory as long as the CEF excitation energy is much 
higher than thermal energy. In fact, the CEF excitation 
energies higher than a room temperature were suggested 
for Np02 5,6 and UO2 9, 10 from neutron scattering ex- 
periments. The r§ quartet has a T4„ octupole moment 
independent of the dipole moment as mentioned before, 
and we expect non-Curie- Weiss behavior for Np02- On 
the other hand, in the subspace, the octupole mo- 
ment is proportional to the dipole moment [see Fig. 2(b)] . 
Thus, we expect simple Curie- Weiss behavior for U02- 

As shown in Fig. 3(a), the magnetic susceptibility of 
UO2 follows the Curie- Weiss law very well at T > 50 K, 
at least up to 300 K. The deviation from the Curie- Weiss 
law at T < 50 K very likely originates from fluctuations 
relating to the antiferromagnetic transition at 30.8 K. On 
the other hand, in Np02 as shown in Fig. 3(b), the mag- 
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netic susceptibility deviates from the Curie- Weiss law be- 
low 150 K. This temperature is much higher than the 
transition temperature 25 K, and the deviation is un- 
likely to stem from critical phenomena. In addition, the 
phase transition at 25 K in NpC>2 has been argued to be 
due to ordering of T 5u octupole moments, 13-19 which do 
not influence the magnetic susceptibility in the param- 
agnetic phase. As shown in Fig. 3(b), the experimental 
data are fitted well with cq (5). 

Let us now try to explain the validity of the magnitude 
of parameters obtained here. For this purpose, it is use- 
ful to discuss coupling constants among multipoles evalu- 
ated from the fitting parameters 9, 9i, and 02- After sim- 
ple algebraic calculations, we obtain CXd = —145.1 K, 
CA = -38.3 K, and C|A do | = 44.2 K. It is found that the 
order of magnitude of 100 K of the dipole-dipole interac- 
tion is consistent with the value of the spin-lattice relax- 
ation time T\ deduced from 17 0-NMR measurements for 
Np02- 2 ° On the other hand, concerning the magnitudes 
of octupole-octupole and dipole-octupole interactions, it 
is difficult to deduce these values from the experimental 
results. However, for a system with active T 4u octupole 
moments, in general, A and Ado can take significant val- 
ues. In fact, from a microscopic /-electron model on the 
basis of a j-j coupling scheme, 19 ' 21 we have found that 
the magnitudes of A and Ado are comparable to that of 
Ad- Thus, it seems reasonable to obtain the above values 
for the coupling constants among multipoles for Np02- 

It may be claimed that the experimental result for 
NpC>2 can be phenomenologically reproduced by x — 
C I (T — 6) + xo, where xo denotes the so-called van Vleck 
term. However, there seems no reasonable explanation 
why xo is significant for NpC-2 whereas xo is almost zero 
for UO2. In addition, it is unclear how to divide the mag- 
netic susceptibility into "spin" and "orbital" parts in 
/-electron systems with strong spin-orbit coupling. We 
claim that the difference in the temperature dependence 
of the magnetic susceptibility between UO2 and NpC-2 is 
naturally understood by considering the effect of the T 4u 
octupole moments in NpC>2. 

In summary, we have discussed the effects of orbital 
degeneracy on magnetic susceptibility in paramagnetic 
phases. Under certain CEF potentials, we have found 
that orbital or multipole moments coupled to the mag- 
netic moment can modify the Curie-Weiss law. Other 
magnetic properties such as nuclear magnetic resonance 
spectroscopy will also be influenced by orbital or mul- 
tipole moments. For instance, on a simple cubic lattice 
with nearest-neighbor interactions, the ^-component of 
the octupole moment with T 5u symmetry couples to J z 
at a wave- vector q with q x ^ ±q y , n and the magnetic 
response at such wave- vectors will be influenced by the 



r$ u octupole moment even in the paramagnetic phase. 
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